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Roles of turbulence in the context of magnetic reconnection are investigated with special emphasis on the 
mutual interaction between flow (large-scale inhomogeneous structure) and turbulence. In order to evaluate 
the effective transport due to turbulence, in addition to the intensity information of turbulence represented by 
the turbulent energy, the structure information represented by pseudoscalar statistical quantities (helicities) is 
important. On the basis of the evolution equation, mechanisms that provide turbulence with cross helicity are 
presented. Magnetic-flux freezing in highly turbulent media is considered with special emphasis on the spatial 
distribution of the turbulent cross helicity. The cross-helicity effects in the context of magnetic reconnection 
are also investigated. It is shown that the large-scale flow and magnetic-field configurations favorable for the 
cross-helicity generation is compatible with the fast reconnection. In this sense, turbulence and large-scale 
structures promote magnetic reconnection mediated by the turbulent cross helicity. 
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I. INTRODUCTION 

Magnetic reconnection is one of the most important 
processes in plasma physics. It appears in a variety 
of phenomena in geo/astrophysics and fusion physics, 
which include dynamo, solar and stellar flares, magne- 
tospheric substorms, disruptions in fusion devices, etc. 
In the framework of magnetohydrodynamics (MHD), the 
diffusion of magnetic field in magnetic reconnection can 
be described with a basic model proposed by Sweelii and 
Parker,- where conservations of the magnetic fiux, mass, 
and momentum are considered. However, Sweet-Parker 
mechanism of reconnection is too slow to explain typical 
time scales of eruptive processes observed in laboratory 
and in space such as solar fiares. In order to alleviate this 
discrepancy due to the slowness of reconnection, several 
theoretical attempts for the fast reconnection, including 
the Petschek's model,— have been proposed. Turbulence 
is one of the candidate mechanisms that enhance the re- 
connection rate very much. 

Turbulence in magnetic reconnection has been investi- 
gated by several researchers. Analysis of magnetic-field 
observations in the Earth's magnetotail suggested the re- 
connection process in the current sheet has properties of 
turbulence excited by tearing vortices^ The effects of tur- 
bulence on magnetic reconnection were first investigated 
by Matthaeus and Lamkin using numerical simulations of 
pinch sheet. They showed that the electric fluctuations 
combined with the flow and current filament enhance the 
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dissipation in the reconnection zone, and leads to the 
rapid reconnection.^'^ In order to dissolve the discrep- 
ancy between the typical length scale of phenomena and 
the dimension required for the current sheet, a fractal 
nature of the current sheet was suggested.'' Lazarian and 
Vishniac examined the effects of stochastic field in recon- 
nection and concluded that the presence of the stochastic 
components of the magnetic field enhances the reconnec- 
tion rate drastically. They stressed that turbulence alters 
the relevant scales for reconnection almost independently 
of the microscopic dissipation mechanismSiS 

One of the primary effects of turbulence is enhanc- 
ing the mixing or transport of the system, which is rep- 
resented by the notion of turbulent transport such as 
the eddy viscosity, eddy diffusivity, anomalous resistiv- 
ity, etc. For example, in the mean momentum equa- 
tion in a turbulent flow, the eddy viscosity i^t enhances 
the effect of molecular viscosity v as, v ^ v + v-y with 
spatiotemporal variation of v^. The enhanced transport 
due to turbulence leads to destruction of structure or 
reduction of spatiotemporal inhomogeneity. In the real- 
world phenomena, however, roles of turbulence are not so 
obvious or straightforward. Under some circumstances, 
a large-scale inhomogeneous structure is formed or sus- 
tained by turbulence itself. Representative situations are 
seen in turbulent dynamos, where large-scale magnetic 
fields are generated by some statistical properties of tur- 
bulent motioufS^— This can be regarded as an example 
of general situations where breakage of symmetry leads 
to the formation or sustainment of large-scale inhomoge- 
neous structure in turbulence. In the presence of symme- 
try breakage in turbulence, we may expect some effects 
that balance or counteract the primary turbulence ef- 
fects. Here, enhanced transports due to turbulence are 
suppressed, and large-scale inhomogeneous structures are 
sustained or formed. Important question is how to for- 
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mulate such breakage of symmetry in turbulence which 
leads to the structure formation. 

Most straightforward method for describing turbulent 
motions is to directly solve governing equations. How- 
ever, for realistic turbulence of scientific interests with a 
huge Reynolds number, it is impossible in the foreseeable 
future to numerically solve all scales of motions rang- 
ing from the largest scales of system size to the smallest 
scales where the energy is dissipated into heat. In such a 
situation, modeling of turbulent motions provides a use- 
ful tool for investigating realistic turbulence. In practice 
of studying turbulence phenomena that appear in the 
geo/astrophysical and fusion sciences, one of the impor- 
tant tasks is how to choose quantities that appropriately 
describe the statistical properties of turbulence. 

In the studies of turbulence modeling, among sev- 
eral statistical quantities, the turbulent energy density 
(per unit mass) k{= (u'^)/2) and its dissipation rate 
e[= i'{{du'°' /dx^)"^)] have often been adopted as descrip- 
tors of turbulence properties (u': velocity fluctuation, 
v. kinematic viscosity).!^ The turbulent energy k ex- 
presses how much fluctuation a turbulent flow has, and 
the energy dissipation rate e expresses how much energy 
is transferred from large to small scales. In this sense, 
both k and e represent the intensity property of turbu- 
lence. The eddy-viscosity representation is supplemented 
by a model of whose generic form is expressed by in- 
tensity measures of turbulence as — /{fc,e}. The 
most popular expression is = C^k^/e {Ci,: model 
constant, k/e: time scale of turbulence). However, it 
is well known in the engineering field that such a simple 
eddy-viscosity model completely fails if applied to some 
turbulent flows. The turbulent swirling flow is the rep- 
resentative case. A simple eddy-viscosity model is too 
dissipative to sustain the centerline dent proflle of ax- 
ial mean velocity that is experimentally observed in the 
turbulent swirling flowJ^iii In order to alleviate this dis- 
crepancy, anomalously small value of Ci, has often been 
adopted in the numerical simulations.^^ This fact implies 
that the description based on k and e is not enough and 
we have to take into account another turbulent statistical 
quantity that properly describes the structure property 
or breakage of mirrorsymmetry in the turbulent swirling 
pipe flow.-'^?. 

Large-scale rotation and magnetic fleld are often es- 
sential ingredients in the geo/astrophysical and fusion 
plasma phenomena. If we have large-scale rotational mo- 
tion or magnetic fleld, symmetry of turbulence between 
the directions parallel and antiparallel to the rotation 
axis or magnetic fleld is broken. Such breakage of mir- 
rorsymmetry is expected to be measured by pseudoscalar 
turbulent quantities, which change their sign under the 
inversion of the coordinate system. Since pseudoscalar 
vanishes in a mirrorsymmetric system, a flnite value of 
pseudoscalar indicates the breakage of mirrorsymmetry. 

The turbulent cross helicity is a pseudoscalar which 
represents breakage of mirrorsymmetry, more precisely, 
asymmetry between the directions parallel and antipar- 



allel to the large-scale magnetic fleld. In the presence 
of breakage of mirrorsymmetry in turbulence, the pri- 
mary transports due to turbulence may be suppressed 
by the cross-helicity effects, and large-scale structures are 
formed and sustained even in the presence of a strong tur- 
bulence. These two effects: (i) generation of large-scale 
structures; and (ii) suppression of turbulent transport; 
are two sides of the same coin, the cross-helicity effect. 
Unlike the kinetic and magnetic helicities, the cross he- 
licity can exist even in two dimensional conflguration of 
magnetohydrodynamic plasmas. For physical interpreta- 
tion of the cross-helicity effect and its application to the 
dynamo problems, the reader is referred to Refs. [iTl - fTgl 

Behavior of turbulent statistical quantities relevant 
to the transport coefflcients provides essential infor- 
mation on the evolution of the large-scale structures. 
At the same time, how turbulent statistical quantities 
are spatially distributed and temporally evolve is de- 
termined by large-scale flow and magnetic-fleld struc- 
tures or mean-field inhomogeneities. As is well known, 
the mean-velocity shear coupled with the Reynolds 
stress,— {u'°'u''')dU°' / dx'' , is the main cause for the tur- 
bulence production in a hydrodynamic (HD) flow. In the 
MHD turbulent flow, as will be shown in Sec. [Til we have 
several mean-field inhomogeneities other than the mean 
velocity shear, which may contribute to the energy and 
cross-helicity generations. 

Large-scale structures, which can be denoted in the 
broad sense as flow including the mean flow and 
global magnetic fleld, etc., are determined by turbulence 
through the turbulent transport coeflacients. On the 
other hand, properties of turbulence are determined by 
the large-scale flow structures such as the mean veloc- 
ity shear, vorticity, magnetic-fleld shear, electric current, 
etc. This mutual relationship between the flow and tur- 
bulence is schematically depicted in Fig. [1] 
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FIG. 1. Flow-turbulence relationship. Mean-field structures 
are determined by turbulence through the turbulent trans- 
port coefficients. Turbulence properties are determined by 
mean-field structures through the production rates of turbu- 
lent statistical quantities. 

If we see the magnetic reconnection phenomena from 
this perspective, important points are as follows. Mag- 
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netic reconnection is considered to be promoted by the 
enhanced magnetic difFusivity due to turbulence. If the 
spatiotemporal distributions of the turbulent statistical 
quantities (turbulent energy, cross helicity, etc.) rele- 
vant to the turbulent transports are favorable for the 
magnetic reconnection, the rate of reconnection would 
be increased. The spatiotemporal distributions of these 
turbulent quantities obey the balance of the production, 
dissipation, and transport rates of these statistical quan- 
tities. The production mechanism, the most important 
ingredient, is determined by the large-scale flow struc- 
tures such as the inflow profiles, global magnetic-field 
configurations, etc. If the large-scale structures sustained 
and formed during the magnetic reconnection process are 
consistent with the favorable distributions of the trans- 
port coefficients for the magnetic reconnection, the re- 
connection would be promoted. Otherwise, the recon- 
nection rate should be reduced. Therefore, it is the pro- 
duction mechanism of the turbulent statistical quantities 
due to large-scale structures that should be examined in- 
tensively. In this paper, we focus our attention on the 
turbulent cross-helicity effects, and address the magnetic 
reconnection from the viewpoint of the mutual relation- 
ship between the flow structure and turbulence mediated 
by the cross helicity. 

This paper is organized as follows. In Sec.|TTl the evolu- 
tion equation of the turbulent cross helicity is presented. 
On the basis of this equation, main mechanisms that pro- 
vide turbulence with the cross helicity are examined. In 
Sec. mil the cross-helicity effects in the mean magnetic in- 
duction equation are treated with special emphasis on the 
conditions for the magnetic-flux freezing and the mag- 
netic reconnections. In Sec. IIVI the cross-helicity effects 
in the mean momentum equation are considered. The 
mean flow configuration affects the reconnection rates 
very much. By considering the flow induced by the cross- 
helicity effects, we discuss the influence of the cross he- 
licity in turbulence on the magnetic reconnections. Sum- 
mary and concluding remarks are presented in Sec. El 



II. CROSS-HELICITY EVOLUTION 

A. Reynolds stress and turbulent electromotive force 

Effects of turbulence on the mean flelds in magnetohy- 
drodynamics are represented by the Reynolds (and tur- 
bulent Maxwell) stress 7?. = (u'u' — b'b') in the mean 
velocity (U) equation and by the turbulent electromo- 
tive force Em = (u' x b') in the mean magnetic-field (B) 
equation (b': magnetic-field fluctuation). Note that, in 
this paper, we adopt the Alfven speed unit for the mag- 
netic field etc. (b — b^/y^/TPi b*: magnetic field mea- 
sured in the original physical unit, /i: magnetic perme- 
ability, p: density). In the case of compressible flows, 
we have several turbulent transport mechanisms that are 
intrinsically connected to the density variation. For ex- 
ample, density variation along the mean magnetic field 



gives an important contribution to the turbulent cross- 
helicity generation. Effects of compressiblity or variable 
density can be incorporated into the cross-helicity ar- 
guments in the context of turbulent dynamos and tur- 
bulent transport suppression. However, such arguments 
will make the description much more complicated. At 
low plasma beta, incompressibility is considered to be a 
good approximation although the plasma beta becomes 
about one at the current-sheet region. Considering that 
the essential physics of magnetic reconnection is not so 
different in incompressible and compressible cases, here 
in this paper, as most work following the original Sweet- 
Parker model, we confine our arguments to the incom- 
pressible framework. This treatment never denies the 
possible importance of compressibility in the magnetic 
reconnection study especially in the astrophysical con- 
text. For some examples of the theoretical extensions to 
the compressible or variable-density case, the reader is 
referred to Refs. 20. and |2T1. 

As we see in Appendix A, from a statistical analytical 
theory for inhomogeneous turbulence applied to magne- 
tohydrodynamic (MHD) turbulence, TZ and Em are ex- 
pressed asii 

= (u'^'u'^ - b'^b'^) (la) 
^^Kj^S^'P^mS^'^ + i^MM''^, (lb) 

Em = (u' X b') (2a) 
= -/3J + aB + 7f2, (2b) 

where Kb.{= (u'^ - b'2)/2) is the turbulent MHD resid- 
ual energy, S and Ai are the strain rates of the mean 
velocity and magnetic field, respectively, J the mean 
electric-current density, r2(= V x U) the mean vortic- 
ity. Here, i^k, i^m, ct, /3 and 7 are the transport co- 
efficients whose spatial distributions are determined by 
the properties of turbulence. In Eqs. (|Tbl) and (|2bl) . the 
eddy viscosity and the turbulent magnetic diffusivity, 
and /3, are expressed in terms of the turbulent energy 
K{= (u'2 -f b'2)/2). The Ai-related transport coeffi- 
cient I'M in Eq. (jlbp and the Jl-related coefficient 7 in 
Eq. (|2b[) are expressed in terms of the turbulent cross he- 
licity (velocity-magnetic-field correlation in turbulence) 
W{= (u' • b')). Then, it is important to examine the 
evolution equations of the turbulent energy and cross he- 
licity, K and W. In the context of turbulent transport 
suppression and structure formation, in particular, the 
production mechanisms of W should be intensively in- 
vestigated. 

Turbulent transports such as the eddy viscosity, turbu- 
lent resistivity, dynamo coefficients, etc. are determined 
by statistical properties of turbulence. A simplest way 
to treat such statistical properties is adopting one-point 
turbulent statistical quantities and consider the evolution 
of these statistical quantities. In magnetohydrodynamic 
(MHD) turbulence, the turbulent MHD energy K and 
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the turbulent cross helicity W defined by 



W : 



(3) 



(4) 



are most important ones among several statistical quan- 
tities since they are directly connected to the transport 
coefficients appearing in TZ, and vm, and the ones in 
Em, P and 7. For detailed description of the transport 
equations of several statistical quantities, the reader is 
referred to Refs. ^ and[23l. 



B. Equations of turbulent MHD energy and cross helicity 

In order to see how and how much energy and cross he- 
licity exist in turbulence, we have to consider the trans- 
port equations for K and W. From equations of fluctu- 
ation velocity and magnetic fields, u' and b', equations 
for the turbulent MHD energy (density) and turbulent 
cross helicity (density) are written as 



sg+Tg 



(5) 



with G = {K,W). Here Pg, sg, and Tg are the pro- 
duction, dissipation, and transport rates of G. They are 
defined by 



Pk = -n''' 



du'^ du"' 
dx°- dx°- 



.dip 
dx°' 



Em • J, 



db"' db"' 
dx°' dx°- 



= B • WW V • T' 



if' 



w 



ah 



dx"- 



Em • S^, 



Tw = B-WK + W -T'n 



(6a) 



(6b) 



(6c) 



(7a) 



(7b) 



(7c) 



(77: magnetic diffusivity) . In Eqs. ([Sa|) and ([7a|) . TZ and 
Em are the the Reynolds stress and the turbulent electro- 
motive force defined by Eqs. ([Ta|) and (l2al) . respectively. 
In Eqs. ([5c)) and (fTc)) . and are the higher order 
correlations of fluctuations whose exact expressions are 
suppressed here. 



C. Cross-helicity generation 

We see from Eq. ([5]) with Eq. ([7]) that there are three 
main mechanisms that provide turbulence with the cross 
helicity. They are the effects related to (i) the Reynold 
stress, —'R!^^[dB^/dx°')] (ii) the turbulent electromotive 
force, —Em • ^\ and (iii) inhomogeneity along the mean 
magnetic field, B-Vif. Among these three, the Reynolds- 
stress- and turbulent-electromotive-force-related produc- 
tions are associated with the drain of the mean-field cross 
helicity U • B. In this sense, they are related to the cas- 
cading property of the cross helicity. On the other hand, 
the turbulent cross-helicity generation due to the inho- 
mogeneity along the mean magnetic field is not directly 
related to such cascade processes but arises from a kind 
of transport with asymmetry of the boundary conditions. 
In this sense, the inhomogeneity along the mean magnetic 
field may play a special role in the turbulent cross-helicity 
generation in the context of externally injecting cross he- 
licity. 

In the following, we refer to some typical situations in 
which each of these generation mechanisms may work. 



1. Production due to Reynolds stress 

The Reynolds stress coupled with the mean mag- 
netic shear gives a cross-helicity generation. Substituting 
Eq. (Hal) into the first term of Eq. (TTaj) . we have a pro- 
duction due to the Reynolds stress as 



b\2 



(8) 



Since the sign of vu is the same as that of W as we see 
from Eqs. (jA8l) and (jAllI) . the second or fM-related term 
always destroys the cross helicity regardless of the sign 
of W. So, as the production mechanism of the turbulent 
cross helicity, we write only the i^j^-related term as 



OB 



1 



(9a;° 2 



(9) 



This shows that a positive cross helicity is generated 
if the mean velocity and magnetic shears have the same 
signs. A typical situation may be caused by a neutral 
beam injection (NBI) in the toroidal direction to toka- 
mak devices in fusion plasmas. In toroidal confinement 
plasmas, the toroidal magnetic field decreases as the ma- 
jor radius of plasma increases (i?*^ oc R: distance 
from the major axis). If a neutral beam is injected in 
the center (minor axis) region, the toroidal velocity is 
maximum at the center of the minor radius (r: distance 
from the minor axis) and decreases as the minor radius of 
plasma increases. In this situation, we have negative ve- 
locity and magnetic shears {dW^/dr < 0, dB'I'/dr < 0) in 
the half region far from the major axis (far side), and neg- 
ative velocity and positive magnetic shears {dW^/dr < 
0, dB'^/dr > 0) in the half region near the major axis 
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(near side) as shown in Fig. [2l This suggests that we 
have positive and negative cross-helicity generations in 
the outer- and inner-half regions of plasma, respectively. 



NBI 




poloidal 



FIG. 2. Cross-helicity production due to the neutral beam 
injection (NBI). The Reynolds stress coupled with the mean 
magnetic-field shear leads to the positive or negative cross- 
helicity generation. 



2. Production due to turbulent electromotive force, 

The turbulent electromotive force coupled with the 
mean vorticity H gives a cross-helicity generation. Sub- 
stituting Eq. (|2b[) into the second term of Eq. ((7a)) . we 
have a production due to the turbulent electromotive 
force as 



Em • = /3J • O - -aB -n. 



(10) 



Since the sign of 7 is the same as that of W as seen in 
Eq. (jAlip , the second or 7-related term in Eq. (|10p works 
for the destruction of the turbulent cross helicity regard- 
less of the sign oi W. In tokamak plasmas, we have a 
very strong plasma current Jq- In such a case, the tur- 
bulent cross-helicity generation due to the electromotive 
force may be expressed as 



(11) 



since |Jo| ^ \^\- This shows that the positive or negative 
turbulent cross helicity is generated if the mean vorticity 
is aligned with the mean electric current: 



Jo • rj > Pw >0, 
Jo • rj < ^ Pw <0. 



(12) 



In tokamak plasmas, we have a strong plasma current 
in the toroidal direction, (2: toroidal direction in 
the cylindrical approximation for a torus geometry). In 
the presence of poloidal rotation , we have a mean 
toroidal vorticity associated with this poloidal rotation, 
fi^[= {l/r){d/dr)rU'^]. Then, it is highly probable that 
we have a generation of positive or negative turbulent 
cross helicity depending on the direction of poloidal ve- 
locity (Fig.El). 

As for the cause of poloidal rotation, we may consider 
several spontaneous plasma rotation associated with the 
improved confinement mode, such as local poloidal ro- 
tation near the plasma edge in H mode, global poloidal 





FIG. 3. Cross-helicity production due to the coupling of the 
large-scale electric current and vorticity. If the large-scale 
vorticity associated with the poloidal rotation is parallel (or 
antiparallel) to the plasma current, a positive (or negative) 
cross helicity is generated. 



rotation near the internal transport barrier in the nega- 
tive magnetic shear mode, etc. If we use a NBI with a 
tilted angle, we have a poloidal rotation in addition to the 
toroidal rotation. However, direct observation of the tur- 
bulent cross helicity in high temperature fusion devices is 
still very difficult. To validate this type of cross-helicity 
generation, a numerical simulation of cylindrical geom- 
etry mimicking the tokamak plasmas was performed;^ 
In the work, it was shown that a poloidal plasma ro- 
tation, which corresponds to the mean vorticity in the 
toroidal direction, coupled with the plasma current con- 
tributes to the cross-helicity generation in turbulence. It 
was also shown that if the direction of poloidal rotation is 
changed, the sign of the generated cross helicity changes. 



3. Production due to inhomogeneity along the mean 
magnetic field 

If the turbulent MHD energy K{= (u'2 + b'2)/2) is 
inhomogeneously distributed along the mean magnetic 
field B, the cross helicity may be generated in turbulence. 
A positive (or negative) cross helicity is produced for V K 
parallel (or antiparallel) to B: 



B VK >0 Pw >0, 
B S/K <0 Pw <0. 



(13) 



It is ubiquitous in astrophysical phenomena that an am- 
bient magnetic field threads through a plasma gas disk. If 
intensity of turbulence is not homogeneously distributed 
in the direction perpendicular to the disk, we have in- 
homogeneity of turbulence along the mean magnetic 
field. This mechanism may play a very important role 
in cross-helicity generation in astrophysical turbulence.?^ 
We consider a situation where a large-scale magnetic field 
threads through a disk from above to below, and the in- 
tensity of turbulence is highest at the midplane and de- 
creases with the distance from the midplane (Fig. 2]). In 
such a case, we have positive and negative cross-helicity 
generations in the upper- and lower-half regions of the 
disk. Consequently, we have positive and negative distri- 
butions of the cross helicity in the upper and lower halves 
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of the disk as 



W>0 
w <0 



for z > 0, 
for z < 0. 



(14) 




p^>0 W>0 
Pw<0 W<0 



FIG. 4. Cross-helicity production due to the energy inhomo- 
geneity along the large-scale magnetic field. 

We should stress that, in this situation, the intensity 
of turbulence is maximum at the midplane so that the 
gradient of turbulent energy, V K, is zero and changes its 
sign there. As this result, there is no cross-helicity gen- 
eration at the midplane, and the magnitude of the cross 
helicity should vanish there. This is consistent with a 
basic feature of pseudoscalars, representing the breakage 
of mirrorsymmetry in turbulence. 



III. CROSS-HELICITY EFFECTS IN MAGNETIC 
INDUCTION 

A. Magnetic-flux freezing in turbulence 

The primary effect of turbulence in the magnetic induc- 
tion equation is destruction of magnetic field, which is ex- 
pressed by the turbulent magnetic diffusivity or anoma- 
lous resistivity /3. If we have no field-generation effects 
such as the helicity or cross-helicity effect, the turbulent 
electromotive force is expressed as 



Em = (u' X b') = -/3J. 



(15) 



Substituting this into the mean induction equation, we 
have 

— =Vx (UxB)-Vx [(r; + /3)VxB]. (16) 

This shows that the effective magnetic diffusivity is en- 
hanced by turbulence as 



T] ^ J] + /3. 



(17) 



We should note that at large magnetic Reynolds num- 
ber (Rm S> 1), the turbulent magnetic diffusivity is in 
general much larger than the molecular counterpart as 



P/r] Rm > 1. 



(18) 



In the absence of any field-generation mechanisms, the 
mean induction equation is written as 



dB 
'dt 



In this case, the magnetic-flux freezing to the plasma 
motion does not occur since f3 is so large that the effective 
magnetic diffusivity can not be neglected. 

On the other hand, if we have some field-generation 
mechanisms such as the helicity and the cross-helicity 
effects, the turbulent electromotive force is expressed as 

Em = (u' X b') = -/3J + aB + 7J7. (20) 

Then, the mean induction equation reads 

— V X (U X B) - V X [(77 + /?) V X B] 
at 

+Vx(aB + 7J7). (21) 

The third term on the right-hand side of Eq. (pij) rep- 
resents the field-generation mechanisms due to the tur- 
bulent residual helicity coupled with the mean magnetic 
field, aB, and the turbulent cross helicity coupled with 
the mean vorticity, jft. If such field-generation mecha- 
nisms are large enough, they may balance the turbulent 
magnetic diffusivity effect, /3J, as 



aB -|-7f2 - 133. 



(22) 



In this case, the contribution from the turbulent electro- 
motive force is very small, and we have the mean induc- 
tion equation as 



(}B 

V X (U X B) + r/V^B. 



(23) 



V X (U X B) - V X [/3V X B] , 



(19) 



This equation has exactly the same form as the induction 
equation for the instantaneous magnetic field b. Then, 
for the flow with very high magnetic Reynolds number 
{Rm ^ 1), the mean magnetic field B is subject to the 
first term of the right-hand side of Eq. 1^^. Again, we 
have a magnetic-fiux freezing to the plasma motion. 

Although this situation looks very similar to the 
magnetic-fiux freezing in the laminar case, it is com- 
pletely different from the laminar one. First, in Eq. ()23|) . 
we have very strong turbulence or fluctuating compo- 
nents of the velocity and magnetic flelds. This magnetic- 
flux freezing occurs only at very high Reynolds and mag- 
netic Reynolds numbers flow. In this sense, we may call 
this the turbulent magnetic-flux freezing. Second point 
is that since this state is based on an equilibrium condi- 
tion (|22p . the "turbulent magnetic-flux freezing" highly 
depends on the mean-field configurations and spatial and 
temporal distribution of turbulent quantities. If the bal- 
ance relation between the field destruction and genera- 
tion breaks, the picture of magnetic-fiux freezing can not 
hold any more. 

If the field-generation mechanisms are absent, the field- 
destruction due to turbulence prevails. A typical situa- 
tion may occur on the midplane of a plasma gas disk or 
equatorial plane in a spherical convection zone. Recalling 
the fact that pseudoscalar turbulent quantities are de- 
scriptors of breakage of mirrorsymmetry, it is quite natu- 
ral these pseudoscalar quantities are distributed antisym- 
metricaly with respect to the midplane. For such spatial 
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distributions, the magnitudes of pseudoscalar quantities 
vanish on the midplane, and very small in the vicin- 
ity of it. As a consequence of this disappearance of 
psedoscalars, a strong effective magnetic diffusivity may 
cause the magnetic reconnection in the close vicinity of 
the midplane. 



nient for both the magnetic-flux freezing in some region 
and the magnetic reconnection in other region. In this 
sense, it is quite natural that model simulations that uti- 
lize the ft effect are very sensitive to the variation of 
velocity profile including the incorporation of the merid- 
ional rotation and the value and spatial distribution of 
the effective magnetic diffusivity. 



B. SI and cross-helicity effects 
1. n effect 

In the study of dynamo, the so-called ft effect has been 
considered as a basic ingredient of the dynamo process. 
The O effect contains several different physical processes 
in it. In the context of the toroidal-field generation from 
the poloidal magnetic field, at least three processes: (i) 
magnetic-flux freezing; (ii) appropriate differential rota- 
tion; (iii) magnetic-field reconnection; should be included 
in the effect. 

As has been already mentioned in Sec. IIIIl in order for 
a large-scale magnetic-flux freezing to occur in a highly 
turbulent medium, we need mechanisms that balance or 
cancel the strong turbulent magnetic diffusivity effect. 
Without such balancing or transport suppression mech- 
anisms, we only have Eq. (|19p. then no magnetic-flux 
freezing occurs. The cross-helicity effect, as well as the 
helicity or a effect, is one of such balancing mechanisms. 

In addition to the magnetic-flux freezing, in order to 
get a toroidal component of the magnetic field from the 
poloidal magnetic field, we need a differential rotation 
of the toroidal velocity. The toroidal velocity should be 
faster or slower in the lower-latitude region than in the 
higher-latitude region. Without such differential rotation 
profiles of the toroidal velocity, the toroidal magnetic 
field can not be generated from the poloidal magnetic 
field in the framework of the 51 effect. The toroidal ro- 
tation observed at the surface of the Sun, the prograde 
rotation in the equatorial or lower-latitude region, is con- 
sidered to support the scenario of a dynamo that utilizes 
the fl effect. 

During the process generating the toroidal magnetic 
field from the poloidal magnetic field by differential ro- 
tation, we in general need a magnetic-field reconnection, 
which makes a new configuration possible from the old 
one. Original poloidal field is distorted by the differ- 
ential rotation. Since such distortion is most strong at 
the midplane or equatorial plane, it is most likely that 
the reconnection occurs there. Then we have a toroidal 
magnetic-field configuration whose direction is opposite 
with respect to the midplane. 

These considerations immediately suggest that the 51 
effect utilized for the toroidal magnetic-field generation 
from the poloidal one is based on a subtle combination of 
several requirements. It requires a particular type of dif- 
ferential rotation profile that is suitable for the toroidal 
field generation. It requires a particular spatial distri- 
bution of the effective magnetic diffusivity that is conve- 



2. Cross-helicity effect 

If the cross-helicity effect is the main balancer for the 
turbulent magnetic diffusivity effect, the mean induction 
equation in a stationary state has a special solution as 
follows. 

We divide the mean magnetic field and electric-current 
density as 



B = Bo + (5B, J = Jo + <5J, 



(24) 



where Bq and Jo(= V x Bq) are the mean magnetic field 
and electric-current density without the cross-helicity ef- 
fect, while (5B and 53{= V x (5B) are the mean fields rep- 
resenting the first-order contribution through the cross- 
hehcity effect. Substituting Eq. (HH) into Eq. ^ with 
the a-related term dropped, we have equations 



at 



V X (U X Bo) - V X (;9V X Bo) (25) 



for the zeroth-order field, and 

V X (U X (5B) - V X (/3V X (5B - 7V X U) (26) 



at 



for the first-order field. 

Equation (j26p represents the magnetic induction 6Ji 
arising from the turbulent cross helicity. Equation (j26p 
has a particular solution 



T W 



(27) 



for a stationary state. Here, Cw(= C^/Cfi) is a model 
constant whose value is estimated as O(10~^) — 0(1) [See 
Eqs. (|XTO)) and (|XTT|) ]. Here, the proportional coefficient 
W/ K is the turbulent cross helicity scaled by the turbu- 
lent MHD energy. This scaled cross helicity is the most 
important quantity in the cross-helicity dynamo formu- 
lation. The value of W/K is not positive-definite, but its 
absolute value is bounded as 



M 

K 



l(u'.b') 



u'- 



b'2)/2 



< 1. 



(28) 



The cross-helicity dynamo solution [Eq. (P7|] requires a 
finite value of turbulent cross helicity for it to work. In 
the absence of cross helicity in turbulence, Eq. (|27p gives 
no contribution at all for the magnetic-field generation. 
In this sense, the spatial and temporal distributions of 
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the cross helicity should be intensively examined. In or- 
der to fully investigate such distributions, we have to 
solve the evolution equation of W [Eq. ([5|) with Eq. d?])] . 
However, without considering the full complications of 
the cross-helicity evolution, we can derive several con- 
sequences of Eq. ([27| as follows. The turbulent cross 
helicity is a pseudoscalar, which describes breakage of 
mirrorsymmetry. As has been referred to in Sec. I, cross 
helicity is considered to be a measure of the asymme- 
try with respect to the magnetic field. If a magnetic 
field threads through plasmas with the turbulent energy 
distributed symmetrically with respect to the midplane, 
we have asymmetry between the directions parallel and 
antiparallel to the magnetic field. Then, we see from 
Eq. (fT4|) that the spatial distribution of the turbulent 
cross helicity is antisymmetric. If such an antisymmetric 
distribution of W is coupled with the symmetric toroidal 
velocity, from Eq. ()27|) we have antisymmetric toroidal 
magnetic field (Fig. [5]). 



Equatorial 
plane 




W<0 



FIG. 5. Toroidal magnetic-field generation due to the cross- 
helicity dynamo. An antisymmetric toroidal-field configura- 
tion is induced from a symmetric velocity profile coupled with 
an antisymmetric cross-hehcity distribution. No differential 
rotation is required. 



A poloidal magnetic field coupled with the inhomoge- 
neous turbulence leads to antisymmetric distribution of 
the cross helicity in turbulence [Eqs. ([T3l) and (|14|) ]. This 
turbulent cross helicity coupled with the toroidal velocity 
gives the toroidal magnetic field [Eq. ((27)) ]. A prominent 
feature of this scenario is that it does not require a par- 
ticular differential rotation of the toroidal velocity for the 
sake of toroidal magnetic-field generation. In the frame- 
work of cross-helicity dynamo, differential rotation is not 
directly related to the toroidal field generation. Differ- 
ential rotation is important only to sustain the energy 
and cross helicity in turbulence. As we see in Eq. p8|. 
the magnitude of the turbulent cross helicity is bounded 
by that of the turbulent MHD energy. Unless we have 
enough amount of turbulent MHD energy, we will not 
have enough cross helicity in turbulence. Only in this 
sense, the differential rotation of plasmas is required for 
the dynamo process to work. 

Very recent numerical simulations showed that a dy- 
namo model constituted of the evolution equations of the 
toroidal and poloidal magnetic fields as well as the turbu- 
lent cross helicity can reproduce an oscillatory butterfiy 
diagram of the solar magnetic activity without resort to 
the n effect 1^ These results will be reported in the forth- 
coming paper. 



IV. CROSS-HELICITY EFFECTS IN MAGNETIC 
RECONNECTION 

Magnetic reconnections are associated with large-scale 
plasma flows with some typical configurations. We focus 
our attention on such large-scale flow structures associ- 
ating with the magnetic reconnection, and examine how 
turbulence influences such flow structures as well as how 
such flow structures determine the properties of turbu- 
lence. 



A. Cross-helicity effects in momentum balance 

So far we treated only the mean induction equation. 
Here we examine the cross-helicity effects in the equa- 
tion of motions. For this purpose, we consider the mean 
vorticity equation 



on 

'dt 



= V X 



V X 



u 



1b 



X n -I- i-kV^ 



u-2b 



1 1 
-(UxB)xB---xB 



(29) 



which is a direct consequence of the expressions for the 
Reynolds stress TL [Eq. ([Tb]) ] and the turbulent electro- 
motive force Em [Eq. ([2b|) ] substituted into the mean 
velocity equation (F: mean external force. A: mean vec- 
tor potential). For the derivation of Eq. (|29)) . see Ap- 
pendix |b] 

As we have seen in Ec^. ([5]) with Eq. ([7]), it is difficult to 
sustain a finite cross helicity in turbulence in the absence 
of the mean magnetic field. In order to clearly under- 
stand the role of the turbulent cross helicity in momen- 
tum transport, we divide the mean velocity and vorticity 
as 



u = Uo + (5u, n = cto + m. 



(30) 



Here, Uq and JIq are the mean velocity and vorticity 
without the mean magnetic-field effects, while 5\5 and 
5Vt are the mean fields representing the first-order effects 
of the mean magnetic field through the turbulent cross 
helicity. Substituting Eq. (|30]) into Eq. (l29|), we have 
equations 



= V X (Uo X J^o + i^cV^Uo -f-F) (31) 



for the zcroth-order field, and 



dsn 



SV 



B X flo + i^kV^ 5V 



/3 



B 



(32) 

for the first-order field. 

Equation ([3T|) expresses the evolution of the large-scale 
fiow structure rig, which is subject to the turbulent vis- 
cosity t'K and the external force F. On the other hand, 
Eq. ([5^ represents the variation of the large-scale flow 
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structure Sfl arising from the turbulent cross helicity or 
7 effect. We see from Eq. ([5^ that, for a given large-scale 
flow structure flo, 



7 W 
p K 



(33) 



is a particular solution of the stationary state. Equa- 
tion p3p suggests that we have a variation of momentum 
transfer that is proportional to the large-scale magnetic 
field multiplied by the normalized cross helicity. 

Here we should make some remarks on the validity of 
Eq. ([55)1 in the compressible case. In compressible flows, 
we have several other terms arising from the density vari- 
ation in Eq. p9| or (jB7p . If we neglect the density fluctu- 
ation (p' = 0) and consider only the mean density strati- 
fication p = p (p: mean density, p': fluctuation density), 
one of the main contributions is the baroclinic effect in 
the mean vorticity equation: 



V X 



4vp 

p 



Vp X VP 



(34) 



If we further assume that the cross helicity does not affect 
the mean density and pressure (p = pg and P = Pq or 
Sp = SP = 0), Eq. dSH) is reduced to 



V X 



P J 



Vpo X VPo 



Po 



(35) 



and no contribution to the first-order equation [Eq. (j32p ]. 
In such a case, the special solution [Eq. ([33]) ] is obtainable 
even in the compressible case. 



B. Reconnection environment viewed from the large-scale 
flow structure 

It is pointed out that properties of magnetic reconnec- 
tion such as the scaling law and maximum reconnection 
rate depend heavily on the form, nature, and value of the 
imposed boundary conditions j^i^ Here we examine the 
reconnection environment from the viewpoint of large- 
scale structure. We consider the circumstances or global 
surroundings of the reconnection region are highly turbu- 
lent. Such turbulence is a consequence of the large-scale 
inhomogeneity of the mean fields such as the velocity 
strain, electric current, etc. In this sense, turbulence 
in magnetic reconnections here is driven by large-scale 
structures. At the same time, we also consider the self- 
generation mechanisms of turbulence in the magnetic- 
reconnection process. 

Since the argument leading to the result [Eq. (l33l) ] is 
quite general, we can consider any general situations in- 
cluding the magnetic reconnection in three dimensions. 
For the sake of simplicity of description, here we consider 
a simple situation where a current sheet is associated with 
a large-scale inflow as in Fig. Ufa). 

Recalling the cross-helicity generation mechanism dis- 
cussed in Sec. Ill Cl we see that the generation mechanism 
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FIG. 6. Magnetic reconnection and turbulent cross helicity. 
(a) Magnetic-field lines and streamlines for a large-scale in- 
flow conflguration. Solid lines: magnetic field. Dashed lines: 
streamline. Shaded region represents the reconnection or dif- 
fusion region, (b) Directions of the mean fields and induced 
velocity due to the turbulent cross helicity. Directions of the 
mean electric-current density J and the mean vorticity Jlo are 
shown with the sign of the turbulent cross helicity W . Direc- 
tions of the induced velocity due to the cross-helicity effect 
[Eq. (|33|l ] are also indicated by the dashed arrows. 



related to the turbulent electromotive force [Eq. (|10p ] 
may be most important in this flow configuration. As 
we see in Fig. [BJb) , the large-scale vorticity J7o is spa- 
tially distributed antisymmetrically with respect to the 
midplane of current sheet (x axis, plane with y = 0), 
and also with respect to the center plane (y axis, plane 
with a; = 0). Since the large-scale electric-current density 
J is aligned with the negative z direction everywhere in 
this case, the sign of the product of Oq and J, J • rJoi is 
distributed as 



J • J7o < for xy> 0, 
J • Jlo > for xy < 0. 



(36) 



Consequently, the production of the turbulent cross he- 
licity is positive and negative in the regions with < 
and xy > 0, respectively. Then it is highly probable that 
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the turbulent cross helicity is spatially distributed as 



Consequently, the magnetic field 



14^ < for a;?/ > 0, 
W > for xy <0 



(37) 



B = Bn + 



(40) 



as depicted in Fig. IHl^b). 

Outside of the current sheet region, it follows from 
Eq. (p3| that the velocity variation due to the cross- 
helicity effect is expressed as 



(5U = 



\w\ 

-Cw-rr-B for x > 0, 
K 



Mb = Cw^B for x<0, 
p K 



(38) 



From the symmetry of the plasma flow, the large-scale 
vorticity vanishes on the midplane (y — 0) and center 
plane (x = 0), so we have no cross-helicity- generation 
mechanism there. The turbulent cross helicity changes its 
sign across the current sheet (midplane) and the center 
plane. For instance, in the positive x region in Fig. |51 
the sign of W changes from positive {y > 0) to negative 
{y < 0) across the current sheet. As this consequence, the 
magnitude of the cross helicity is small in the vicinity of 
the current sheet. Equation ([38|) thus makes a very small 
contribution to the velocity in the current-sheet region 
(and also in the center- plane region). 

In the absence of the turbulent cross helicity {W — 
7 = 0), we have no induced flows, then the flow conflg- 
uration remains as it was originally set (Uo and i~lo). If 
the flow is a potential one, we have no mean vorticity. In 
this situation, we have no cross-helicity generation since 
J • flo = everywhere. A stagnation-point-like flow is 
such an example. On the other hand, in the presence of 
the turbulent cross helicity, the induced flow due to the 
cross-helicity effect should alter the flow proflle. How 
much flow configurations are modulated depends on 7//? 
or W/K, the turbulent cross helicity scaled by the turbu- 
lent MHD energy. For example, in the case of maximum 
cross helicity (|VF|/i^ = 1), the induced flow speed is 
close to the Alfven speed. As this result, outside the 
current-sheet region, the inflow is converging toward the 
central plane (x — 0), and in the current-sheet region the 
outflow is concentrated very narrow region in the vicinity 
of the midplane (y = 0), where the magnitude of W is 
very small. 

In the mean magnetic equation, the turbulent cross 
helicity coupled with the large-scale vortical motion gives 
rise to the induced magnetic fleld that is aligned with the 
mean velocity [Eq. (P7|) ]. The spatial distribution of the 
turbulent cross helicity, Eq. ([57)) . leads to the induced 
field: 



SB 



_Mu = -Cw^U for xy > 0, 
p K 

Mu = Cw^U forxy<0. 



(39) 



turns to be a curved field configuration as schematically 
depicted in Fig. [T] Combination of the inflow converg- 
ing toward the y axis with the curved magnetic fleld 
line gives a X-point-like configuration, which is known 
to be suitable for the fast reconnection. In the absence 
of the turbulent cross helicity, such modulations of the 
large-scale velocity and magnetic fields never occur. This 
suggests that the presence of cross helicity with a favor- 
able mean-field configuration leads to the enhancement 
of magnetic reconnection rate as compared with the case 
without cross helicity. 



W>0 



W<0 




W<0 



W>0 



FIG. 7. Flow and magnetic-field configurations in fast mag- 
netic reconnection. Solid lines: magnetic field, Dashed lines: 
plasma fiow. Shaded region represents the reconnection or 
diffusion region. Signs of the turbulent cross helicity W asso- 
ciated with the slow shocks are also indicated. Sign of W is 
negative (or positive) for the region where the Alfven wave is 
propagating in the directions parallel (or antiparallel) to the 
large-scale magnetic field. 



C. Reconnection rate 

As we have already mentioned, the width of the diffu- 
sion region in magnetic reconnection is reduced because 
of the converging flow due to Eq. ([55)1 and the curvature 
of the magnetic field due to Eq. ([39t . which would lead 
to the enhancement of the reconnection rate. Here, we 
estimate the reconnection rate with the aid of the expres- 
sions for the cross-helicity effects in the magnetic induc- 
tion [Eq. (|27p ] and in the momentum balance [Eq. ([55]) ]. 

As the dimensionless reconnection rate, we adopt the 
inflow Alfven Mach number— 



Mi, 



Uin 



(41) 



where is the magnetic field measured in the Alfven 
speed unit (Bin = Vmu = BlJ ^fjip, VAin: Alfven speed, 
B;*-,: magnetic field measured in the physical unit). 

Let us denote the inflow velocity and magnetic fleld in 
the absence of the cross-helicity effects as ul"' and B^^' , 
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respectively. From Eqs. (1331) and ([27)) . the inflow velocity 
and magnetic field with the cross-hclicity effects, Uin and 
Bin, are expressed as 



(0) j_ 7 



Bin = b1"' + SB., = B 



(0) ^ 7 



Ui, 



(43) 



respectively (Fig. |5]). This system of equations are writ- 
ten in a matrix form as 



1 -A^u. 



(44) 



where *Uin = (C/f„,C/f„), *Bin = [Bf^^Bl), etc. Solving 
Eq. p4)) with respect to *Uin and *Bin, we have 



*Ui, 
*Bi, 



H -1 



1 2\ /tTT(0)\ 



where m/j^-* = U^^' /Bia is the reconnection rate in the 
absence of the cross-helicity effect, and 6*0 is the angle 



AO) 



between u|'°'' and Bj^'' . In the case of converging inflow 
and X-point like curved magnetic field such as in Fig. [7j 

(42) 

the sign of cos^q or IJ^^^ ■ Bj^'' obeys the quadrupole 



,(0) 



distribution similar to that of W [Eq. ([57)1 ] 

{cos ^0 < for xy > 0, 
cos ^0 > for xy < 0. 



(49) 



Then, the sign of (27//3)Mi^°^ cos 6*0 in Eq. (gS]), as well 
as the (7//3)^-related contributions, is always positive. In 
this sense, in the converging inflow cases it is not the 7//? 
and cos Oq themselves but the magnitudes of them which 
determine the modulation of Min. 

If U|°^ is perpendicular to b|°' {Oo = 7r/2), Eq. dH]) is 
reduced to 



MM = ^/2) 



M 



(0)^ 



(7//?)' 



(7//3)^m(°) 



1 



1/2 



(50) 



Then, for given U|°^ and b|°\ we can calculate Uin and 

Bin. 



dUi, = (y//5)Bi, 




= iyip) Ui, 



The behaviors of Afin with respect to |7|//3 with several 
values of in the case ofM/°) =0.1 are plotted in Fig. [HI 
Here, as an example, we consider the domain with xy > 
0, where 7 < 0. As expected, we have no enhancement of 
Min by the cross-helicity effect if 7//3 = 0. With increase 
in I7I//3, Min increases up to 1 (at |7|//3 = 1). Here, we 
should note that 7//? is expressed by the turbulent cross 
helicity scaled by the turbulent MHD energy as 



FIG. 8. Modulation of the velocity (Left) and magnetic field 
(Right) due to the cross-helicity effects. In this figure, the 
case with the negative cross helicity {W < 0, 7 < 0) is drawn. 



From Eq. (j45|) . the magnitude of the inflow velocity 
and magnetic field are given as 



Uin = lUinl 



-1/2 



u 



(0)^ 



m 111 



1/2 



(46) 



-1/2 



in in 



J J m ■ m ■ ^ 

Then, the reconnection rate is expressed as 



1/2 



M„ 



M^f + (7//3)^ + C^l/mC cos^o 
(l/f^yM^f + 1 + (27//3)m(°) cos 00 



(47) 



1/2 

, (48) 



M 

/3 



— Cw 



m 

K 



< 1 



(51) 



with 1 > Cw = 0(10 ^)— 0(1), which corresponds to the 
fact that [Eq. (UTTj) ] is smaller than Cp [Eq. ([MO]) ] 
but the orders of them are the same. 

These results suggest that the magnetic reconnection 
rate may be substantially enhanced by the cross-helicity 
effects if the magnitude of the scaled turbulent cross he- 
licity is large enough such as 



m 

K 



> 0.1. 



(52) 



On the other hand, the cross-helicity effect may be neg- 
ligible for the magnetic reconnection rate in the case of 
the turbulent cross helicity of 



M 

K 



< 0.1. 



(53) 



In the case of non-converging inflow geometry, the 
signs of W and cos may be opposite to each other, 
then the sign of {2j//3)M^-^^ cos 6*0 is negative. In such a 
case, the cross-helicity effects work for reducing the re- 
connection rate for small |7|//3-value domains since the 
contribution of (27/;3)m/°^ cos ^0 is more dominant than 
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FIG. 9. Variation of tlie reconnection rate Min against tiie 
scaled turbulent cross helicity 7//?. The angle between the 
inflow velocity and magnetic field in the absence of the cross- 

helicity effects, UjjJ' and b||J', is denoted by ^0. : ^0 ~ 

7r/2; : So = 27r/3; - - -: So = 57r/6; , So = tt/S; 

- •• -, So = 7r/6. In all plots, the magnetic reconnection rate 
without the cross-helicity effects, M^^\ is set equal to 0.1. 



ity rJo associated with the plasma flow is distributed as 
depicted in Fig.[TUl Then the production of the turbulent 
cross helicity due to the coupling of J and flo is expressed 
as in Eq. ([55)1 , leading to the turbulent cross helicity spa- 
tially distributed as in Eq. ([37]) . These points should be 
carefully scrutinized using numerical simulations. 



y 




FIG. 10. Tearing-mode configuration and distribution of the 
turbulent cross helicity. 



V. CONCLUDING REMARKS 



that of (7//?)^. In Fig. [HI the Mjn behavior in such situ- 
ations {60 = 7r/6, 7r/3) are also plotted. 

Another important point to remark is that the turbu- 
lent cross helicity distributions preferable to the recon- 
nection are expected to be sustained in the down-stream 
region or at the after-reconnection stage. If there is an 
asymmetry between the directions parallel and antipar- 
allel to the magnetic field, we have a certain production 
mechanism of the non- vanishing turbulent cross helicity. 
As is well known, X-point reconnection and slow shock 
regions are source of the Alfven waves. For instance, 
the anisotropic ion beams in the plasma sheet boundary 
layer, which is the so-called PSBL ions in the magneto- 
tail,^* is known to excite the Alfven waves by the ion cy- 
clotron beam instability^^ The reflected ion beams from 
the slow shock front may excite the Alfven wave in the 
shock upstream along the magnetic-field line as welli^ 
Since the Alfven waves generated by slow shocks propa- 
gate in the outward direction, the sign of the cross helicity 
associated with the Alfven waves is negative (or positive) 
if the large-scale magnetic field is parallel (or antiparal- 
lel) to the direction of the wave propagation. In terms of 
turbulent statistical quantities, this mechanism is related 
to the cross-helicity generation due to the inhomogeneity 
along the mean magnetic field, B • WK [Eq. ([T^]. The 
resultant cross-helicity distribution is consistent with the 
one expressed by Eq. ([57)) as shown in Fig. [71 

In relation to the spatial distribution of the turbu- 
lent cross helicity, we should note that the quadrupole 
symmetric distribution of the turbulent cross helicity 
[Eq. (j37p ] may be commonly or ubiquitously observed in 
the magnetic-reconnection configurations. For example, 
in the tearing-mode configuration, the large-scale vortic- 



Magnetic reconnection was viewed from the flow- 
turbulence interaction mediated by the cross-helicity ef- 
fects. In addition to the enhancement of transport, the 
primary effect of turbulence, turbulent transport may be 
effectively suppressed by some other effects of turbulence. 
The structure information of turbulence represented by 
some pseudoscalar statistical quantities such as the tur- 
bulent cross, kinetic, current, magnetic helicities is the 
key for analyzing the transport suppression or structure 
formation due to turbulence, whereas the intensity in- 
formation of turbulence is essential for describing the 
transport enhancement or structure destruction. In or- 
der to investigate the balance between these two effects, 
we flrst considered the cross-helicity production mecha- 
nisms in turbulence. The cross-helicity effects, which ap- 
pear in the Reynolds stress [Eq. ([Ibp ] and the turbulent 
electromotive force [Eq. (|2bp ]. were examined in the con- 
text of the magnetic-flux freezing in turbulence, turbu- 
lent dynamo process, and magnetic reconnection. In the 
arguments on magnetic reconnections, we stressed the 
importance of the large-scale flow structure and global 
magnetic-field configuration since these inhomogeneous 
structures constitute the conditions for the cross-helicity 
generation. The spatiotemporal distribution of the cross 
helicity, in turn, determines the transport properties of 
turbulence, then affects the magnetic reconnection rates. 
We showed that the large-scale flow and magnetic-fleld 
structures associated with the magnetic reconnection are 
favorable for the cross-helicity generation in some cases, 
and that the flow and magnetic-field structures induced 
by the cross-helicity effects are consistently favorable for 
the enhanced magnetic reconnections. In the absence of 
the turbulent cross helicity, these mutual interactions be- 
tween the large-scale structure and turbulence will never 
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appear at all in the magnetic reconnection process. 

This theoretical scenario of magnetic reconnection 
based on the mutual interaction between the large-scale 
structures and the turbulent cross helicity should be val- 
idated using numerical simulations. A few points should 
be remarked on such simulations. First, the theoretical 
formulation based on the cross-helicity effects through 
Eqs. (|27p and ([55)) extended in this paper is very general, 
and not limited to the two-dimensional cases. Estimate 
of reconnection-rate modulation and numerical valida- 
tion of it in three dimensions are desired much especially 
in the context of astrophysical magnetic reconnection. 
However, we can start with two-dimensional cases since, 
unlike the kinetic and magnetic helicities, the cross helic- 
ity can exist even in the two dimensional cases. Secondly, 
boundary conditions for the numerical simulations may 
be important but in a limited sense. Actually, it was 
reported that the conditions on flows are important in 
determining several features of magnetic reconnection, 
including the nature of the inflow, the length of the dif- 
fusion region, etci^ However, in the scenario presented 
in the present paper, such boundary conditions are im- 
portant only in the sense that they affect the spatial dis- 
tribution of the turbulent cross helicity (more generally, 
distributions of relevant turbulent statistical quantities) 
through which properties of turbulent transport are de- 
termined. Numerical simulations validating the present 
theoretical idea are being prepared, and results will be 
reported in the forthcoming papers. 
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Appendix A: Cross-helicity dynamo and its solution 

1. Turbulent electromotive force 

From a closure theory for inhomogeneous turbulence 
applied to magnetohydrodynamic (MHD) turbulence, 
the Reynolds (and turbulent Maxwell) stress Ti- and the 
turbulent electromotive force Em are known to be ex- 
pressed asii 

7^"^ - ^i^R - vkS''^ + v^M'^P, (Al) 



where Kb.{= (u'^ — b'^)/2) is the residual energy, and 
<S and Al are the strain rates of the mean velocity and 
magnetic fields defined by 



M 



al3 _ 



(A3) 



(A4) 



In Eqs. (|A1|) and (IA2p . i^k, i^m, P, 7, and a are the 
transport coefhcients, which are known to be expressed 
as 

a ^ - J dk J dTiG{k,x;T,Ti,t) X 

X [-Hy^uik, x; T, Ti,t) + Hbb{k, x; r, Ti,t)] , (A5) 



f3 — - J dk J dTiG{k, x; r, ti, t) x 

X [Q„„(fc,x;T, Ti,t) + Qb6(fc,x;T,Ti,t)] , (A6) 



j^-jdkj dTiG{k,x]T,Ti,t) X 

X [(3„b(fc,x;T, Ti,t) + Qbu{k,x;T,Ti,t)] , (A7) 



z.K = (7/5)/3, I'M -(7/5)7. 



(A8) 



Here, G is the Green's function of MHD turbulence, and 
Hum Hbb, Quu, Qub, Qub are the spectral functions of the 
turbulent kinetic helicity, current helicity, kinetic energy, 
magnetic energy, and cross helicity. Equations (jA5[) - (jA7[) 
show that a, /3, and 7 are determined by the spectral 
distributions of the residual helicity, energy, and cross 
helicity, respectively. These coefhcients have spatiotem- 
porally non-local nature through the Green's function. 

The simplest possible models for the transport coef- 
ficients a, /3, and 7 are given by using turbulent one- 
point quantities: the turbulent residual helicity H{= 
(-u' • b' • j')), the turbulent MHD energy K{= 
(u'^ + b'^)/2), and the turbulent cross helicity W{= 
(u' • b')) respectively: 



7 = C^tW 



(A9) 



(AlO) 



(All) 



(Ca, Cp, and C^: model constants). Here t is the time 
scale of turbulence, which is often expressed by using the 
turbulent MHD energy scaled by its dissipation rate: 



E 



M 



-/3J + 7^ + aB, 



(A2) 



K/e. 



(A12) 
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2. Cross-helicity dynamo solution 

If we drop the a- or helicity-related term in the turbu- 
lent electromotive force Em [Eq. (jA2p ] as 



Em - -l3J + -fn, 
the mean induction equation (I2ip yields 
dB 



(A13) 



dt 



Vx(U X B)+Vx(-/3V X B + 7V X U), (A14) 



where the molecular magnetic diffusivity rj have been 
dropped as compared with the turbulent counterpart /3. 
In a stationary state, the mean magnetic field 



7 W 



(A15) 



is a special solution of Eq. (|A14I) provided that the spatial 
variation of /3/j{— CwW/K) is not so large.— 

Appendix B: Vorticity equation in the presence of the 
turbulent cross-helicity effects 

The mean momentum equation in the rotational form 
is written as 
9U 



dt 



= Uxr2 + JxB-V-7?. + F 



-V P 



1, 



1. 



.'2 



(Bl) 



where 7?. is the Reynolds stress defined by Eq. (Hal) . F 
the mean external force, P the mean pressure function. 
If we substitute expression (fTb]) for TZ into Eq. (|B1[) . we 
have 



dV 

'dt 



U X J7 + J X B- 



-V ( P + 



i^KV^ ( U - ^B I + F 



(B2) 



Then, the mean vorticity equation is given by 
U X f2 + J X B 



-i/kV^ U 



1 



B 



- V X F. 



(B3) 



In addition to the Reynolds stress, which contains fiuc- 
tuating velocity correlation and the fluctuating Lorentz 
force (j' X b'), the turbulence effects enter the mean 
velocity equation also through the mean Lorentz force 
J X B. The mean Ohm's law is written as 



J = o-(E + U X B + Em) . 



(B4) 

If we substitute expression (I2b[) for the turbulent electro- 
motive force into Eq. (jB4l) and solve it with respect to J, 
we have the mean electric-current density as 



J==i(UxB 



o 



(B5) 



(A: vector potential). Here, we have dropped the molec- 
ular magnetic diffusivity 77 as compared with the turbu- 
lent one P (rj <^ /3). Then, the mean Lorentz force is 
expressed as 

JxB = i(UxB)xB + ^f2xB-i^xB. (B6) 

Here, we should note that the 7- or cross-helicity-related 
term remains in the mean Lorentz force expression. This 
gives a strong contrast to the a- or helicity-related term, 
which disappears due to the alignment of aB with J in 
Eq. (|B5[) . In other words, unlike the turbulent helicity, 
the turbulent cross helicity contributes to the mean ve- 
locity or vorticity equation through the mean Lorentz 
force. 

Substituting Eq. (jB6[) into the mean Lorentz force in 
Eq. (|B3p . we obtain the equation for the mean vorticity 
as 



dfl 

'dt 



= V X 



V X 



U - ^B ) X n + i^kV' 



u-1b 



1 1 f)\ 

-(UxB)xB---xB 



(B7) 



This mean vorticity equation was originally derived in 
the context of system rotation effect As long as the 
expansion (jSOp is adopted and the induced vorticity is 
considered, this equation can be also utilized for cases 
without a system rotation. 
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